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Abstract 

In 1993, a proof was published, within this journal, that there are no regular solutions to 
the linearized version of the twisting, type-N, vacuum solutions of the Einstein field equations. 



While this proof is certainly correct, we show that the conclusions drawn from that fact were 
unwarranted, namely that this irregularity caused such solutions not to be able to truly describe 
pure gravitational waves. In this article, we resolve the paradox — since such first-order solutions 
must always have singular lines in space for all sufficiently large values of T — by showing that if 
we perturbatively iterate the solution up to the third order in small quantities, there are accept- 
able regular solutions. That these solutions become flat before they become non-twisting tells 
us something interesting concerning the general behavior of solutions describing gravitational 
C^) , radiation from a bounded source. 
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1. Introduction. 

The generic behavior of gravitational radiation from a bounded source is clearly an 
important physical problem. Even reasonably far from that source, however, type-N solu- 



^ ■ tions of the vacuum field equations must have both non-zero twist as well as reasonable 



asymptotic behavior in order to provide an exact description of that radiation. Such so- 
lutions would provide small laboratories to better understand the complete nature of the 
singularities of type-N solutions, and could also be used to check numerical solutions that 
include gravitational radiation. It is therefore reasonable that there is considerable interest 
in this problem. In addition to non-zero values for the twist parameter, interesting solu- 
tions must also have appropriate asymptotic behavior; the only currently-known solution 
to the twisting problem, due to Hauser [1], does not have this asymptotic behavior. As 
well the definitive relevance of homothetic vectors to the study of metrics of Petrov type 



N has been very well enunciated by Mcintosh [2]. It is the existence of two homothetic 



vectors — more exactly, an H 2 of symmetries, a non-Abelian group of (local) nomothetic 
vectors for the manifold — that allows the defining equations to be reduced to an ordinary 
differential equation. By construction, this equation must contain a constant parameter, 
the homothetic parameter, which Mcintosh called n. Hauser's solution has the value of 
5/2 for the Mcintosh parameter, n. 

Although several distinct formulations of this problem already exist, none have yet 
been able to produce new solutions. We therefore believe that any exact solutions to this 
problem would be important in obtaining a better understanding of the general problem. 
Clearly Hans Stephani felt much the same way when he published his proof [3] that all, 
non-flat, first- order solutions of the twisting, type-N, vacuum Einstein field equations must 
always contain singular lines, i.e., places on the £, ^-sphere such that for sufficiently large 
values of r, and u, the analytic functions contained within the solution must a) be non- 
constant since otherwise the solution would be flat, and b) therefore must contain poles on 
that sphere which extend to infinite values of the affine parameter, r. Since it is true that 
such singular behavior is the rule for non-twisting Petrov type N solutions, Stephani then 
conjectured that all pure, Petrov type N solutions of the vacuum field equations would be 
insufficient to completely describe the propagation, in vacuum, of gravitational waves from 
a bounded source. We resolve this apparent paradox, in the discussion below, by defining 
an algorithm that allows us to extend Stephani's discussion of the original Einstein field 
equations to an arbitrary order. We then find that at third-order we can indeed show the 
existence of a solution that is non-flat (of Petrov type N), twisting, and also regular for 
sufficiently large spheres. 

2. The description of the vacuum field equations: 

To most easily compare the equations to other forms, we give our presentation in the 
second variant of the usual null tetrad formalism as originated by Debney, Kerr, and Schild 
in Section 4 of their paper [4]; however, for ease of comparison to the work of Stephani, 
we (mostly) use the same symbols as he does, which come originally from the work of 
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Kramer, Stephani, MacCallum and Herlt [5]. Therefore we write the metric, g, in terms 
of a complex null tetrad, as follows: 

g = 9^®e v = 2e 1 ®e 2 + 2e 3 ®e 4 , e 1 = e 2 , e 3 , e 4 real, (2.1) 

s s s 

and use the overbar for complex conjugation. For any vacuum, type-N space-time with 
non-vanishing complex expansion, i.e., where Z = — r 42 i 7^ 0, Debney, Kerr, and Schild 
showed the existence of local coordinates {£, £}, complex, and {r, it}, real, with r the affine 
parameter along the radiation trajectories. In terms of these coordinates, they showed 
that one can always write the null tetrad, e M , so that 

e 1 = ^-dC , e 2 =^=d(, e 3 = du + Ld( + Ld( , 

PZ PZ (2.2) 

e 4 = dr + Wd( + Wd( - He 3 , 

where the metric functions are given by 

Z~ x = (r-iZ), 2zE = P 2 (DL — DL) 

W = -±L iU + iDZ, D = d c -Ld u (2.3) 

H = -rd u (\nP) + \K , K = 2P 2 Re[D(D\nP - L u )] 

where the subscripts indicate partial differentiation. Within this tetrad, setting P = V u , 
the remaining Einstein vacuum field equations take the form 



D^P~ l d u DDV j = , (2.4a) 
DDDDV = DDDDV . (2.46) 

Contrariwise, we insist that the solutions be non-flat, and have non-zero twist, which insists 
that both the following two quantities should not vanish: 



CW oc d u d u {P^DDV} ^ , 

(2.5) 

2i E =P 2 (DL - DL) ^ . 
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3. The description of the perturbation equations: 

To define the perturbation scheme, we first notice [5] that one may always take a 
gauge condition designed to maintain the function P at its simplest value, namely 



p = i + RC • (3.i) 

As well, if we now append the requirement that L = I/(C>C> u ) = 0' we S e t exactly the 
following form of the flat space-time metric: 

ds 2 = 2{ T — — \ d(® dC + 2 du® dr + du® du , (3.2) 

[I+2CCJ 

where ( and ( simply define the (usual) stereograhic coordinates on the complex sphere, 
while r is the radial coordinate and u is the retarded time coordinate. Our perturbative 
scheme can then be developed by thinking of L as being defined in terms of a repeti- 
tive sequence of approximations, involving higher and higher order approximations, and 
defining 



$ = -P~ X DDV = L c + ^-=L-Ld u L , (3.3) 

1 + |CC 



(_ 

_ 1 
2 

in terms of them. Therefore, Eq. (2.4a) now take the form 

Dd u $ = = &^d u <$> - Ldl$ . (3.4) 

It is now obvious that = implies that $^°- ) = 0, so that we may start with the 
lowest-order approximation, denoted in this form. To find the next step — the first-order 
approximation — we suppose that both L and $ may be expanded into series, and take 
Eq. (3.3) as defining $ for us in terms of L, and then take Eq. (3.4) as the constraint 
which defines $ within the next order, so that we immediately have the constraint 

dzd u *W = =► ^ 1 )=«W(C,«) + /5 (1) (C,C) • (3-5) 
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Inserting this expression back into Eq. (3.3), we obtain a defining relation for namely 

d c L™ + — = $U + L^d u LV + L^d u L^ = a^((, u) + /?«(C, () . (3.6) 
1 + 2 CC 

The general solution of this equation, for L^ l \ is given by 

^ (1) = (i + KO- 2 {/ (1) (C,«) + /rfC(i + K0 2 [« (1) (C,«) + /? (1) (C,0]} • (3.7) 

This allows us to write, now, the second-order steps in the form 



d L (2) + _^_ L (2) = $(2) + L W duL W . 

i + Kc 

The general solution of these equations is then easily given in the form 
$(2) = j du j dCLMdl$M , 

= (1 + KO" 2 {/ (2) (C, u) + / d( (1 + KO 2 [* (2) + L^d u L^] } 
At the n-th level, the equations to be solved are then just 



n-l 



L (n) + I <_j L (n) = $ (n) + £ L (n~j) Q^U) . 



n-l 



(3.8) 



(3.9) 



(3.10) 



The general (n-th order) solution of the equations for $ and L are then given by the series, 
up to the n-th order, of the $( J ) and the determined as 

n-l 



*<»>= J du J dCL {n - 3) d 2 u ^K 

L (n) = (1 + 1^-2 | /( n) ( ^ u) + J d( (1 + 1^)2 ^(n) + ^-,^0)] } 



(3.11) 



It now remains to consider Eq. (2.4b) which, in the current notation, can be written 



Im(DDP^) = 
5 



(3.12) 



Expanding out the terms in D, we obtain 



Pirn j<%%$ + - L<9 U <%$ - L (d u $r$ - Ldl§) 

1 + 2"CC 



( chL + %-=L - Ld, 

V 1 + ice 



< u Ljd u $}=0. 



(3.13) 

However, we can use Eq. (3.3) and Eq. (3.4) to re-write substantially the terms in the 
above equation, giving us 

Im + - ^- d-$ - $d u <f> - Ldfdjb | = . (3.14) 

Therefore, applying the same procedures as we have been using so far, for a constraint on 
the n-th iterative step, we find the linear, inhomogeneous, pde for <3>( n ), namely 

Im jd^$ (n) + #c$ (n) j = Yl lm {® {n - j) d u ® U) +i>- J )<%d u $ (i) } • (3.15) 

This equation surely does have solutions, and we may conclude that our iterative procedure 
is indeed complete. Within the n-th order approximation, one has that 

n n 

$(n) ' L ~ Yl Lin) ■ ( 3 - 16 ) 

4. The first few approximate solutions 

Returning, first, to the linearized, or first-order approximation, we of course have 
that the solution is given by Eq. (3.5) and Eq. (3.7), where = a^\(,u) and = 
f^'iCi u )i while the function (3^ = fl^iC, C) is allowed to be an solution of the equation 



Im 



|^a) + _^=^a)| = o. (4.i) 



Defining, now, the function E = E((, u) such that 

d?E((,u)=a( 1 \(,u) , (4.2) 
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we can quickly re- write in the form 



^'^ (l + Ko 2 2(i + |co 2 i + Kc 

where B(C,C) = (j^T^) / rf C (1 + |CC) 2 (C, C) , 



and we have denoted our f^\(,u) by Stephani's symbol C((,u), while we have used 
E(£,u) to denote his -D(C, «)• With these identifications, this is exactly the solution given 
by Stephani in [3]. Insisting that the solution be non-flat, we infer that it is necessary for 
the existence of a non-trivial, linearized Petrov type N solution that we must have 



dldfE^O . (4.4) 

However, as has already been pointed out by Stephani [3], this can happen if and only if 
the field is singular at some point on the £, ^-sphere. Every non-singular field defines a flat 
spacetime, only, at this level of approximation. 

To see this in some more detail, we note, still following Stephani, that the invariants 
available to us are 

p2 

(4.5) 



^-^ (i) =i-{^ (i) -^ (i) }, 

J2 « + J 2 (1) = 1 {l - P 2 (d c d u L {1) + HlW)} . 

Differentiating the first of these with respect to u, and comparing the result with the 
second, one easily concludes that if both of these are regular on the £, ^-sphere, then 
it must also be true that P 2 d^d u L^ must be regular there. Employing our equation, 
Eq. (4.3), for we easily calculate that this quantity is given by 

P%d u L^ = -—^d u C(C,u)+d^d u C(Cu) + —^d u E(Cu)-d < d u E(Cu) . (4.6) 

We easily infer that this quantity is, then, regular on the sphere, if and only if the 
functions C and D are of the form 

C(C,u) =(r(u) + a(u)+p(0, 

(4.7) 

E((,u) =(5(u)+u;(u) + h(0. 
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However, the form given in Eq. (4.7) for £"(£, it) would obviously lead to the space-time 
being flat, rather than, non-trivially, of Petrov type N. This result, given by Stephani [3], 
seems to lead to the following conclusion, which we refer to as Stephani's 'paradox,' namely 
[3] 

So either twisting type N fields do not describe a radiation field outside a 
bounded source, or due to a mechanism not yet recognized, they cannot be lin- 
earized, or the naive interpretation of the coordinates (starting from their New- 
tonian limit) is wrong. 

The most catastrophic of these options is the first possibility, i.e., that "twisting, type N 
fields do not describe a radiation field outside a bounded source." In the next section, we 
find a solution for the twisting, type N field, at the third step of iteration. Our solution 
is in fact regular on the £, ^-sphere. This allows us to resolve the catastrophic argument 
cited above, by saying that it is indeed not true. 

5. Example of a regular solution at the third level of approximation. 

Referring back to Eq. (4.3) for l/ 1 ), we choose it to have the very simple form 

a (l) a (l) ( 7 i~] 
L (1) = — — , = — — < - t , a (1) a complex constant. (5.1) 

It is straightforward to check that this solution does indeed satisfy Eq. (3.15) for n = 1. 
Now, we may choose — see Eq. (3.9) — the second-order quantities to be given by 

Once again, it is evident that the constraint Eq. (3.15) is satisfied for n = 2. Now, using 
the n-th order equation for n = 3, we find that 

L ' 3 ' = (TTW { ?/<3,(u) + I dC t (1 + ^ )24><3> + ~'" m ^1 1 ' (5 ' 3) 



We may then set as yet another complex constant, and put 

* (3, = — t + -{tt^-?}- (b ' 4) 

Once again, we may substitute this value for <E>( 3 ) into Eq. (5.3) for which gives 
us its form: 

"(TTW + 7 ^TTlc? + ^- (5 ' 5) 

Inserting, as before, these values into Eq. (3.15), we find that they satisfy that equation 
for n = 3. 

We must now determine the values for the twist and the curvature that go along with 
this, third-order approximation for the complete solution. We have, immediately, that the 
twist, E, is given by 

£(°) = 0, 

E (D = Ipa (^ L (D _ a^ 1 )) = o, 

E( 2 ) = lp 2 (^)-^) = ^Im(/( 2 )(.)), (56) 

E (3) = Ip^ ^ L (3) _ ^^(2) _ + ^(1)^^) j 

= [^Im(/(3)( w ))- 2 [Re(aW)]Im(& 

along with the iterated values for the contributions to the (2-dimensional) curvature, K: 
K (0) = 2P 2 Re(DDlnP) = 1 , 
= -2P 2 Re(d c d u L {1) ) = 0, 

i + ICC V du J (57) 

= -2P 2 Re(d ( d u L {3) - L^dlL {2) ) = - 2 1 "f C £ Re ( ^ u) ) 

u 1 + |CC V du J 



so that the total, 4-dimensional curvature, CW = 2\I>4, is given iteratively by 

$ (°) = , tf^ = , *f } = , 

, T ,(3)_^ A + KC \ V/ (2) N ( 5 - 8 ) 
4 r V C / ' 

Bringing together the sum of all these quantities, to write the complete solution valid 
up through the third iteration step, we have 

^ _ ai ( (_ A ai df 2 (u) 

~ c ^i + |cc c) e du ' 

^-^(TTW^^TTi^— ■ (59) 



1 + -|CC \ du 

K«l- 2 i4«Re + 4|Re(oi)]Re (*hW 



I+5CC \ du J V 
along with the curvature itself, 

where we have defined 

O! = a« + a< 2 > + a< 3 > , f 2 (u) = f^ 2 \u) + /< 3 >(«) . (5.11) 

In order to actually determine the metric itself, we must, lastly, determine the function, 
W, which, to this level of iteration, is given by 

W r— h jlm/ 2 (u) - — f-= r — — -~ h iha- 



(l + iCC) 2 V du J ^>) Cl + KCV du* du 

(5.12) 

Inserting all this into the equation for the metric, g, itself, we see that the metric does 
indeed appear to be everywhere regular on the £, ^-sphere. To completely show this, one 
needs only to change the coordinates in a neighborhood of the north pole, according to 
the usual rule, (' = and ( = 1/C> which causes no trouble at all. 
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Gathering all our results together, we conclude that indeed one can find a regular, 
twisting, non-flat, Petrov type N, vacuum metric that is regular on the £, C -s phere 
in the third order of approximation. Consequently it seems that the twisting, type N fields 
can describe a radiation field outside a bounded source. Of course it is quite interesting to 
determine what happens in the next iteration steps. We intend to consider this question 
soon. 
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